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Abstract 

This is a follow-up to our earlier work for the energies and the charge (vector) and matter (scalar) 

distributions for S-wave states in a heavy-light meson, where the heavy quark is static and the light 

quark has a mass about that of the strange quark. We now study excited states of these mesons 

with higher angular momentum and with radial nodes. 

The calculation is carried out with dynamical fermions on a 16^ x 32 lattice with a lattice spacing 

I a ~ 0.11 fm. The lattice configurations were generated by the UKQCD Collaboration. Attempts 

psj ' are now being made to understand these results in terms of the Dirac equation. In nature the 
> . 

. closest equivalent of this heavy-light system is the Bg meson, which allows us to compare our 

O 

I lattice calculations to experimental results (where available) or give a prediction where the P-wave 

' states should lie. We pay particular attention to the spin-orbit splitting, to see which one of the 
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(N 



o 



X 



states (for a given angular momentum L) has the lower energy. 
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I. MOTIVATION 



There are several advantages in studying a heavy-light system on a lattice. Our meson 
is much more simple than in true QCD: one of the quarks is static with the light quark 
"orbiting" it. This makes it very beneficial for modelling. On the lattice an abundance of 
data can be produced, and we know which state we are measuring - the physical states 
can be a mixture of two or more configurations, but on the lattice this complication is 
avoided. However, our results on the heavy-light system can still be compared to the Bs 
meson experimental results. 



II. MEASUREMENTS AND LATTICE PARAMETERS 

We have measured both angular and radial excitations of heavy-light mesons, and not 
just their energies but also some radial distributions. Since the heavy quark spin decouples 
from the game we may label the states as L-t = L ± |, where L is the angular momentum 
and ±1 is the spin of the light quark. 

The measurements were done on a 16'^x32 lattice. We have two degenerate quark fiavours 
with a mass that is close to the strange quark mass. The lattice configurations were generated 
by the UKQCD Collaboration. More details of the different lattices used in this study can be 
found in Refs. and references therein. Two different levels of fuzzing (2 and 8 iterations 
of conventional fuzzing) were used in the spatial directions to permit the extraction of the 
excited states. 



III. 2-POINT CORRELATION FUNCTION 

The 2-point correlation function (see Figure [1]) is defined as 

C2{T) = {PtrG,{^,t + T,t)Pt+Tr^U'^{^,t,t + T)) , (1) 

where [/"^(x, + T) is the heavy (infinite mass)-quark propagator and Gg{yi,t + T,t) the 
light anti-quark propagator. Pf is a linear combination of products of gauge links at time t 
along paths P and F defines the spin structure of the operator. The (...) means the average 
over the whole lattice. The energies (m^) and amplitudes (oj) are extracted by fitting the 



2 



nL±-lS 


Lattice 1 Lattice 2 Lattice 3 


nL±-lS 


Lattice 1 


Lattice 2 


Lattice 3 


IS 











2S 


0.68(3) 


0.858(15) 


0.82(3) 


IP- 


0.46(4) 


0.42(2) 


0.43(3) 


2P- 


1.40(6) 


1.38(2) 


1.36(5) 


1P+ 


0.61(4) 


0.50(3) 


0.42(4) 


2P+ 


1.60(5) 


1.44(3) 


1.35(6) 


1D± 


- 


L01(4) 


1.03(6) 


2D± 


- 


1.54(2) 


1.67(5) 


ID- 


0.93(6) 


0.86(3) 


0.86(5) 


2D- 


1.83(4) 


1.771(15) 


1.79(4) 


1D+ 


0.96(8) 


0.91(3) 


0.96(3) 


2D+ 


1.80(6) 


1.80(2) 


1.85(2) 


1F± 


L36(5) 


1.24(4) 


1.20(3) 


2F± 


2.10(2) 


2.029(15) 


2.027(13) 



TABLE I: Heavy-light meson energy differences in GeV. Lattice 1 uses a static lieavy quark, wliereas 
Lattices 2 and 3 are slightly smeared in the time direction. Lattice 2 uses APE type smearing and 



Lattice 3 uses hypercubic blocking. See 



J] for details, a = 0.110(6) fm was used to convert the 



energies to physical units. (Note: In the first version of this article a = 0.115(3) fm was used.) 



C2 with a sum of exponentials, 



max 



C2(T) ^ where iVmax = 2 - 4, T < 14. 



(2) 



i=l 



Fuzzing indices have been omitted for clarity. 



IV. ENERGY SPECTRUM AND SPIN-ORBIT SPLITTINGS 

The energy spectrum obtained is shown in Fig. [2] and in Table [B There seems to be some 
spread in the extracted energy for the P+ state between the different lattices, but the reason 
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FIG. 1: On the left: Two-point correlation function. On the right: Three-point correlation function. 
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FIG. 2: Energy spectrum of the heavy-light meson. Here L+(— ) means that the hght quark spin 
couples to angular momentum L giving the total j = L ± 1/2. 2S is the first radially excited L= 
state. The D+— is a mixture of the D— and D+ states, and likewise for the F-\ — . Energies are 
given with respect to the S-wave ground state (IS). Here a = 0.110(6) fm was used to convert the 
energies to physical units. The error bars shown here contain statistical errors from the lattice 
calculations and the uncertainty in the extracted lattice spacing. The solid line on the one-body 
Dirac equation - see section IVTl for details. 

for this is not understood yet. The energy of the DH — state was expected to be near the 
spin average of the D— and D+ energies, but it turned out to be a poor estimate of this 
average. Thus it is not clear whether or not the FH — energy is near the spin average of the 
two F-wave states, as was hoped. Our earlier results can be found in Ref. jl|. 

One interesting point to note here is that the spin-orbit splitting of the P-wave states 
is small, almost zero. We extracted the energy difference of the P+ and P— states in two 
different ways: 

1. Indirectly by simply calculating the difference using the energies given by the fits in 
Eq. [21 when the P+ and P— data are fitted separately. 

2. Combining the P+ and P— data and fitting the ratio S^tS-^, which enables us to go 

^2(-r -) 

directly for the spin-orbit splitting, ?7ip_|_ — mp_. 

The results of the fits are (—5 ± 45) MeV and (11 ± 13) MeV, respectively. D-wave spin- 
orbit splitting was also extracted in a similar manner, giving the results (102 ±54) MeV and 
(124 ± 18) MeV, respectively. 
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FIG. 3: Interpolation to the b-quark mass. The experimental data is from [3]. For another lattice 
group's nic/mQ = results see [4|. 



To obtain predictions of the meson excited state energies, we can now interpolate 
between the static heavy quark lattice calculations and Ds meson experimental results, i.e. 
interpolate between the static quark and the charm quark. Linear extrapolation seems to 
work well, as can be seen in Fig. [3l The P-wave states seem to lie near the BK thresholds. 



V. RADIAL DISTRIBUTIONS: 3-POINT CORRELATION FUNCTION 

For evaluating the radial distributions of the light quark a 3-point correlation function 
shown in Fig. [1] is needed. It is defined as 

C^{R,T) = {V^U'^VG,,QG,2). (3) 

This is rather similar to the 2-point correlation function in Eq. [1] We have now two light 
quark propagators, Gqi and Gq2, and a probe 6(i?) at distance R from the static quark (74 
for the vector (charge) and 1 for the scalar (matter) distribution). 

Knowing the energies rrii and the amplitudes from the earlier C2 fit, the radial distri- 
butions, x''^{Rys, are then extracted by fitting the C3 with 

^max 

GsiR^T)^ J2 ttie-™'*! x^^(i?) e-™^*2aj. (4) 
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The results are 
errors only. See 



otted in Figs. H] and O The error bars in these figures show statistical 
^ and Q for earlier calculations. 



VI. A MODEL BASED ON THE DIRAC EQUATION 



A simple model based on the Dirac equation is used to try to describe the lattice data. 
Since the mass of the heavy quark is infinite we have essentially a one-body problem. The 
potential in the Dirac equation has a linearly rising scalar part, bscR, as well as a vector 
part byecR- The one gluon exchange potential, «OGE ' ^OGE' modified to 



(5) 



QCD 



where Aqqj^ = 260 MeV and the dynamical gluon mass mg = 290 MeV (see for details). 
The potential also has a scalar term mujL{L + 1), which is needed to increase the energy of 
higher angular momentum states. However, this is only a small contribution (about 30 MeV 
for the F-wave). 

The solid lines in the radial distribution plots are predictions from the Dirac model fit 
with m = 0.088 GeV, aqGE = O-^l. ^sc = 1-14 GeV/fm, foyec ~ f'sc and uj = 0.028. These 
are treated as free parameters with the values obtained by fitting the ground state energies 
of P-, D- and F-wave states and the energy of the first radially excited S-wave state (2S). 
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FIG. 4: On the left: The S-wave ground state charge distribution. The label "model" on the solid 
line refers to the model presented in section IVII On the right: The S-wave ground state and 1st 
excited state charge distribution overlap. Note that we see one node, as expected from the Dirac 
equation. 
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FIG. 5: On the left: The P-wave ground state charge distributions. The P+ distribution has a 
peak a bit further out than the P— , which is expected. The sohd hnes are predictions from the 
model in section IVTl On the right: The D— ground state charge distribution. 



Note that the excited state energies in Fig. [2] were not fitted. 



VII. CONCLUSIONS 

• There is an abundance of lattice data, energies and radial distributions, available for 
this Bg-like system. 

• The P-wave spin-orbit splitting is small (essentially zero), which supports the symme- 
try fovec = ^sc as proposed in [7!]. The D-wave spin-orbit splitting is clearly non-zero 
and positive, whereas another lattice group seems to find it to be slightly negative 



see 



1: 



The energies and radial distributions of S-, P- and D-wave states can be qualitatively 
understood by using a one-body Dirac equation. 
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